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ABSTRACT 

We apply the Bayesian concept of 'evidence' to reveal systematically the nature of 
dark energy from present and future supernova luminosity distance measurements. We 
express the unknown dark energy equation of state w(z) as a low order polynomial 
in redshift and use evidence to find the polynomial order, thereby establishing the 
minimum order required by the data. We apply this method to the current supernova 
data, and with a prior — 1 < w(z) < 1 and VL m — 0.3 ± 0.05, obtain a large probability 
of 91% for the cosmological constant model, with the remaining 9% assigned to the 
two more complex models tested. We also investigate the use of evidence for future 
supernova data sets such as distances obtainable from surveys like the Supernova 
Acceleration Probe (SNAP). Given a low uncertainty on the present day matter density 
we find that, if the underlying dark energy model is only modestly evolving, then a 
constant w(z) fit is sufficient. However, if the evolution of the dark energy equation 
of state to linear order is larger than \w\\ > 0.5, then the evolution can be established 
with statistical significance. For models where we can assume the prior —1 < w(z) < 1, 
the correct polynomial order can be established even for modestly evolving equations 
of state. 

Key words: cosmology:theory - methods: statistical -cosmological parameters. 



1 INTRODUCTION 



Observations of distant supernovae give strong indica- 
tions that the expansion of the universe is accelerating 
JRiess et al. 19981 IFerlmutter et al. 19991 iRiess et al. 200T1 
|Tonry et al. 2003} . Indications for a low matter density 
universe are provided by X-ray observations of rich clusters 
IBahcall fc Fan 19981 |Mohr, Mathiesen fc Evrard 1999| 
IDodelson fc Knox 2000ft . which recently has been confirmed 
with remarkable accuracy ( |Allen, Schmidt fc Fabian 20021 . 
If cosmic microwave background (CMB) observations 
( Hins haw et al. 2003ft are combined with large scale struc- 
ture surveys, such as the two degree Field Galaxy Redshift 
Survey (2dFGRS) (IPercival et al. 20021 . the best fit cos- 
mological model is a flat FRW universe with Qa ~ 0.7 and 
On ~ 0.3 HSpergel et al. 2003ft . 

In the simplest case these observations are explained 
by the addition of a cosmological constant term in the 
Einstein's theory of gravity. The cosmological constant 
is often identified with the energy density of the vac- 
uum (Zel'dovich 1968), but explaining its small value today 
(10 _120 Mp 1 ) in terms of fundamental physics has remained 
unsuccessful l |Weinberg 1989| ICarroll 2001ft . Therefore, at- 
tempts have been made to explain the missing dark energy 



as the energy density of a minimally coupled scalar field 
called Quintessence. 

Due to the slow roll of the Quintessence 
field JWetterich 19881 IPeehles R.atra 19881 

IRatra fc Peebles 19881 the universe can become domi- 
nated by vacuum energy and the expansion begins to 
accelerate. In these models the energy in the scalar field 
becomes important only at relatively late times, giving 
ample time for the growth of structures in the universe; 
while tracker like solutions for the field help towards 
ameliorating the fine tuning problem faced by a pure 
cosmological constant (jStcinhardt, Wang fc Zlatev 19 99 1 . 

At present there are a large number of models that 
can describe the observed acceleration in the expansion of 
the universe. These models differ from a cosmological con- 
stant mainly through their equation of state pq = wpQ. By 
assuming a constant equation of state the current bounds 
are w < -(0.6 - 0.8) jSpergel et al. 2003ft . Although most 
models predict w > —1, this is not necessary in non- 
minimally coupled field models (Uzan 1999 Chib aT999l 
lAmen dola 2000 1 or models with a non-canonical kinetic 
term ( |Armendariz-Picon, Mukhanov fc Steinhardt 2 000 1 . 

One of the most promising probes for dark energy are 
future supernova observations, such as the proposed Su- 
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pernova Acceleration Probe - S NAP 1 ||A ldering et al. 2002| l , 
the w project ESSENCE 2 (IGarnavich et al. 200311 and 
the Canada - France - Hawaii Telescope Legacy Survey 
(CFHTLS) 3 . The potential of a SNAP class experiment in 
distinguishing theoretical models has been investigated in 
great detail in various works. 

Most methods approximate the dark energy equa- 
tion of state w(z) in terms of low order polynomials, 
or other simple fitting functions 
Maor, Brustein & Steinhardt 2001 



(Efstathiou 1999 
IGoliath 20011 



Wcllcr & Albrccht 2002 



Weller fc Albrecht 200lP 

Gerke fe Efstathiou 20021 . A polynomial approxima- 
tion is easily interpreted in terms of weighted means of 
the true equation of state JSaini et al. 20031 . The poly- 
nomial approximation naturally allows us to phrase three 
interesting questions about the dark energy: First is our 
universe described by a cosmological constant (w — —1), 
second is the equation of state a constant and and 
third, the most general case, is w evolving with redshift? 
To answer these questions we need to decide on the basis 
of data which polynomial order is required to fit the data 
adequately, regardless of the fitted values of the parameters. 
The fact that polynomial coefficients are related in a simple 
way to the true equation of state also allows asking further 
questions regarding the possible time evolution of dark 
energy. 

The Bayesian construct 'evidence' answers precisely this 
type of question. The method has the usual likelihood ratio 
test built into it, as well as the sensible criterion for pe- 
nalising models with a large number of parameters, known 
as the Occam's razor term QMacKay 199l1lSivia 19961 . Our 
main aim in this paper is to use this concept to formulate 
the dark energy problem as a model selection problem, and 
to quantify the requirements of the present and the future 
data for elaborate modelling. In addition the method is eas- 
ily extended to differentiate between the diverse theoretical 
models for quintessence, even when their parameterizations 
are totally unrelated. 

In Section [5] we introduce the concept of evidence and 
illustrate its meaning using a simple example. Section[3]con- 
tains a brief overview of Type la supernovae as standard 
candles, and a brief description of the SNAP mission. In 
Section^Jwe use evidence to calculate the acceptance proba- 
bility for simple models of dark energy by using the currently 
available SNe data. We investigate the future prospects for 
discriminating amongst various dark energy models by sim- 
ulating SNAP like data and applying the evidence proce- 
dure. Furthermore we analyse how prior information on the 
matter density and the dark energy parameters is affecting 
evidence. Our conclusions are presented in Section |5] 



2 BAYESIAN MODEL SELECTION: 
EVIDENCE 

The notion of evidence is not frequently used in the astro- 
nomical literature, therefore we review the Bayesian method 



1 See at 

2 See at 

3 See at 
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for model selection, and then apply it to the problem of in- 
ferring the polynomial order required by a given data set. 
For further details see Sivia (1996). 



2.1 Mathematical framework 

A model consists of a set of rules to predict data from a 
given set of parameters and a prior which quantifies the 
probabilities of the different parameter values in the absence 
of any data. Consider a set of models (hypotheses) {H}. 
From Bayes' theorem the probability that H is true is 



P(H\D) = 



P(D | H)P(H) 
P(D) 



(1) 



where D denotes the observed data. This shows how our 
prior probability P(H) is modified by the presence of the 
data to give the posterior probability P(H | D). 

The probability P(D | H) is the probability of data 
marginalized over the parameter values in the model H . This 
can be seen more clearly from the following. Writing the pa- 
rameters for model H as 0, Bayes' theorem gives for the 
posterior probability of the parameters given the data and 
model: 



P(0\D,H) = 



p(p\e,H)P{e\H) 

P(D\H) 



(2) 



P(D\0, H) is the usual likelihood of data, given the model 
and its parameters, and P(6\H) are the priors on the param- 
eters. The required quantity is the denominator in the right 
hand side which is found from normalizing the left hand side 
to unity to be 



(3) 



£ = P(D \H)= d n 0P(D | 6, H)P(0 \ H) 



where £ denotes the evidence of the hypothesis H . 

If the data is predicted by a large volume of the param- 
eter space allowed by the priors then the model gets a high 
probability (evidence) , which is a very desirable feature. The 
denominator in Eqn. Q is an overall normalization constant 
which can be ignored if one is interested only in the rela- 
tive merit of the various hypotheses. The term P{H) is our 
prior probability for the various models being compared. A 
uniform prior over all the models being considered would 
express our lack of inclination for any particular model. In 
this case the posterior probability P(H \ D) in Eqn. for 
the various hypotheses is proportional to the evidence. 



2.2 Interpretation: what evidence measures 

To gain a better intuitive understanding of evidence we now 
carry out the calculations analytically for a simple example. 
Assume that the likelihood P(D | 6, H) is a Gaussian about 
the best fit likelihood position Ol 



P(D\9,H) = P(D \6 L ,H) exp 



- 2 (e-e L ) 



F(6-6 L ) 



(4) 



where F is the usual Fisher (curvature) matrix defined by 
Fij = ~d 2 [log P(D\0,H)] /dOidOj evaluated at 6 = L . As- 
sume that the prior is also a Gaussian, but centred on dp 
with curvature matrix P. Evaluating the integral in Eq. [3] 
we find 
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S = P(D\0 L ,H) exp(-C) 



IF + PI 



-1/2 



(5) 



where C is a constant depending on the degree of overlap 
between the prior and likelihood distributions. For the sim- 
pler case where 0l = Op we find that C = and if we define 
volume under the prior V pr ior = |P| _1/ ' 2 and the volume 
under the posterior V^ ost = |F + P|~ 1//2 then the evidence 
becomes 



Vpost 
Vprior 



(6) 



The first term gives the usual likelihood at the best fit point, 
and the second fraction is the ratio between the volume un- 
der the posterior to the volume under the prior. Therefore 
evidence takes into account the likelihood at the best fit 
point, which usually increases with the number of parame- 
ters. However this is partially countered by the second term 
— the so called Occam's razor term — which ensures that 
models with a large number of parameters are penalised. 



3 SUPERNOVAE 

Type la supernova appear to be an excellent standard can- 
dle (Perlmuttcr ct al. 19971 IRiess et al. 1998(1 . with a small 



dispersion in apparent magnitude, a n 



0.15, with no ev- 



idence for evolution with redshift. The apparent magnitude 
is related to the luminosity distance through 



m(z) = M + 5 log D-l(z) 



(7) 



where M = M + 5 log [(c/H Q ) /Mpc] + 25. The quan- 
tity Mo is the absolute magnitude of Type la SNe and 
Dl(z) = dh(z)/(c/Ho) is the Hubble constant free lumi- 
nosity distance. The combination of the absolute magni- 
tude and the Hubble constant, M, can be calibrated by low 
redshift supernovae, for instance the Calan/Tololo sample 
(Ham uyet al. 1993|IPerlmutter et al. 199911 . The dispersion 
in the magnitude, cr mag , is related to the uncertainty in the 
inferred distance, a by 

a _ In 10 

This is about 7% for (J mag = 0.15. In our simulations we 
assume that the errors in the luminosity distance are Gaus- 
sian. We neglect systematic errors in our calculations. 

The SuperNova Acceleration Probe (SNAP) survey is 
expected to observe about 2000 Type la supernovae up to a 
redshift z ~ 1.7, each year (Aldering et al. 2002). Although 
the expected distribution of SNe is complicated, for our cal- 
culations we assume them to be uniformly distributed. Since 
a single supernova measures the luminosity distance with 
a relative error of ~ 7%, binning the supernovae in ~ 50 
bins would give a relative error in the luminosity distance of 
about ~ 1%. 

For our simulations we calculate the luminosity distance 
in 50 redshift intervals up to a maximum redshift z = 1.7. 
We assume a relative error of 1% in the luminosity distance. 
Note that we do not add noise to the simulated distances. 
Therefore, our results only give the ensemble average of the 
various quantities that we quote below. 

To decide which polynomial order is best suited to the 
data we now apply the notion of evidence. To label different 



models we choose the letter N — the order of the polynomial 
approximation. The case of cosmological constant is denoted 
as N = — 1. The only free parameters are the matter density 
f2 m and the polynomial coefficients w = (wi). In reality we 
also have to marginalize over M , but for simplicity we ignore 
this in our calculations. Therefore, our free parameters are 

= (On,w). 

The likelihood function P(D | 0, N) is given by 



P(D\0,N) = Af f[ exp 



df (zi,G) -&,(*) 



-(9) 



where the index i ranges from 1 to A^at, the number of 
supernovae in our sample (or the number of bins in redshift) 
The data is described by the measured luminosity distance 
dh(zi), the dispersion o~i and the redshift Zj. Also 



dt\z) 



c(l + *) . 
Ho 



(1 + z')- 3 ' 2 dz' 



(10) 



° JaL + n Q exp 1 3 if w(z")/{l + z")dz"} 



where w(z) = ^2f =0 Wizl > ^ m ^ s ^ ne P resen t day fractional 
energy density in pressure-less matter and £Iq = 1 — il m is 
the present day fractional density in dark energy. In Eqn.|5| 
the normalization constant J\f = (27r) _JVdat/2 /n^I d 1 at er>. 



4 EVIDENCE ANALYSIS 

We begin this section by comparing evidence values for dif- 
ferent theoretical models fitted to the current supernova 
data. We then repeat this exercise by simulating future data 
for various assumed input models and show how well evi- 
dence picks out the correct model. We then investigate the 
issues related to the effect of uncertainty in the matter den- 
sity, and the effect of different priors on the dark energy 
parameters. 



4.1 Current supernova data 

Although the quality of data will continue to improve, it 
is instructive to see what the current supernova data tells 
us about which model to trust. For the analysis we use the 
compilation of supernova luminosity distances and redshifts 
by Tonry et al. (2003). We apply a prior of Q m = 0.3 ± 0.05 
(similar to Allen, Schmidt & Fabian 2002). In Fig.Qwe plot 
the 68% and 95% contours. As expected the error bars are 
large. The peak of the likelihood surface lies in the region 
where w < — 1. Table H"Tl gives the computed evidence values 
for various models with two different priors for w(z): the first 
row gives the evidence values with the prior —1.5 < wo < 0, 
—2 < wi < 1 (limits of Figure 0, and the second row gives 
evidence with the prior —1 < w(z) < 1, in the range < 
2 < 1.7. As expected from Figure Q the simplest A model 
has the largest evidence with either prior. The data favours 
the simplest model but still allows the possibility of some 
evolution in agreement with the contours plotted in Fig. 
Since the peak of the posterior does not lie in the region 
given by the second prior — 1 < w(z) < 1, we find that with 
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Figure 1. The 68% and 95% joint likelihood contours for u>o — 
w\ with current Supernova data jTonry et a l. 2003 1, marginalized 
over S7 m with a Gaussian prior, Q m = 0.3 ± 0.05, and a uniform 
prior for the equation of state in the range —1.5 < wo < and 
-2.0 < toi < 10. 



Table 1. Evidence for the Supernova luminosity distances given 
in Tonry et al. (2003). We have used the tight prior Q m = 0.3 ± 
0.05 for this calculation. 



this prior the data is more consistent with the cosmological 
constant. 



4.2 The discriminatory power of future supernova 
surveys 

We simulate luminosity distance for a SNAP like experi- 
ment as described in Section |H] First we consider a linearly 
evolving equation of state given by w(z) = —0.8 + 0.6z 
with f2 m = 0.3. This choice of parameters illustrates cer- 
tain degeneracies first described in Maor et al. (2002) In 
Fig [5] we plot the 68% confidence contours in the wo-wi 
plane, marginalized over f2 m , for different priors on f2 m . For 
the dark energy we chose the priors —1.5 < wo < and 
—2 < wi < 1 in this calculation. We address the issue of 
priors on the dark energy in detail in Section 14.41 We see 
how the joint likelihood contours tighten as Q m is increas- 
ingly constrained. 

We also plot contours using a Gaussian approximation 
to the peak of the posterior for a prior < fi m < 1. Com- 
paring with the corresponding exact calculation shows that 
the Fisher matrix error bars are very misleading. This is not 
surprising since Fisher matrix only gives the minimum vari- 
ance for the extracted parameters. The correct contours also 
show that the distribution of the equation of state parame- 
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-1 



-0.5 



Figure 2. The marginalized 68% joint likelihood contours in the 
Wo — wi plane for an input model w(z) = —0.8 + 0.6z and f2 m = 
0.3 for various priors on Q m - The solid line is for a fiat prior in 
the range < Sl m < 1, the dashed line a Gaussian prior with 
f2m = 0.3 ±0.1 and the dot-dashed line a Gaussian prior with 
f2m = 0.3 ± 0.05. The dotted line is the result of a Fisher matrix 
analysis. 



ters wo and wi is not close to Gaussian. Therefore, for our 
calculations below we do not use Fisher matrix approxima- 
tions described in Sec. 12.21 

We now investigate how well the evidence can infer the 
required polynomial order for the equation of state. For 
this purpose we choose to consider fitting models only up 
to a linear order, including the case of a cosmological con- 
stant model (A model). Therefore, our model space com- 
prises (1) the A model with w(z) = —1; (2) the constant 
equation of state model with w(z) = too and (3) the lin- 
early evolving model with w(z) — wo + w\Z. In Figure |3] 
we plot the normalized probabilities as a function of fit- 
ting models, for various input models. For this figure we 
assumed a tight prior on fi m = 0.3 ± 0.05. Evidence clearly 
picks out the correct polynomial order for the A model (solid 
line), the model with w(z) — —0.7 (short dashed line) the 
model with w(z) = —0.8 + 0.6z (dotted line). Here we have 
also considered a model with a weaker redshift evolution, 
w(z) = —0.8 + 0.3z (long-dashed line). We find that for this 
case the evidence marginally prefers a model with no red- 
shift evolution, although the evidence is very similar for all 
three models. The evolution in the equation of state in this 
model is clearly not sufficient to be established by SNAP. In 
general the evidence tends to be conservative in this respect 
and favours simpler models. 



4.3 Effect of uncertainty in f2 m 

The greatest uncertainty in constraining w(z) comes from 
the fact that the value of fi m is not very well known. To 
illustrate the effect of choosing different priors for f2 m we 
calculate the evidence for model w(z) = —0.8 + 0.6z with 
the following three priors on Q m : (1) A uniform prior in the 
range < fi m < 1, (2) a Gaussian prior with an m ~ 0.1 
and (3) a Gaussian prior with an m ~ 0.05. Our results are 
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Figure 3. Evidence values (in percent) for various models, and 
various levels of fit as indicated on the x-axis, with a Gaussian 
prior f2 m = 0.3 ± 0.05. The solid line is for an input A model, 
the short dashed line for an input model with w(z) = —0.7, the 
long-dashed line for a model with w(z) = —0.8 + 0.3z and the 
dotted line for w(z) = —0.8 + 0.6z. 
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Figure 4. Evidence values for the model ui(z) = —0.8 + 0. 62 with 
various levels of fit as indicated on the x-axis. The solid line is for 
a uniform prior in the range < O m < 1, the dashed line for a 
Gaussian prior, f2 m = 0.3 ±0.1 and the dotted line as in Fig. 151a 
prior of f2 m = 0.3 ± 0.05. 

plotted in Fig. [I] The correct model is picked out only in the 
case where the prior on the matter density is the tightest, 
ffn m = 0.05. 

The real difficulty is the degeneracy that exists between 
a fast evolving model and a cosmological constant model 
with an incorrect value of f2 m IMaor et al. 20021 . This can 
be understood from Fig. |K| where we plot the log-evidence 
for the three models as a function of assumed values of Q m . 
We see that the evidence for the true linear model is nearly 
constant and stays above the evidence for the other two 
models until fi m ~ 0.46, where a degenerate case appears 



which fits best to the cosmological constant model, and the 
worst to the linear model. In fact, the A model with Q m = 
0.46 and the linear model with the correct f2 m have nearly 
the same likelihood, so the evidence decides on the basis 
of the Occam's razor term and penalises the more complex 
model. As a result the uniform prior of < Sl m < 1 prefers 
the A model, and since the Gaussian prior with fi m = 0.3 ± 
0.1 model does not rule out fi m = 0.46 sufficiently strongly, 
it still assigns the largest probability to the A model. In fact 
the correct inference is drawn if either the priors on f2 m are 
tight, or if we have a hard upper bound at f2 m ~ 0.4, as 
Fig. |5] clearly demonstrates. 

This result can also be understood in terms of Fig. 
The cosmological constant model corresponds to the inter- 
section of the axes at wq = —1, wi = 0. For the two widest 
Sl m priors the 68 per cent contours almost include the cos- 
mological constant model, which explains why evidence uses 
Occam's razor to pick this model. Whereas for the tightest 
prior on O m the 68 per cent confidence contours are far from 
the cosmological constant model, and indeed, as expected, 
the integrated probability (evidence) disfavours it. 

Constraints on Q m can be obtained by including other 
cosmological measurements with the supernova data. Some 
probes such as the CMB and cosmic shear (see Refregier 
et al, 2003) are sensitive to a combination of w and fi m , 
so this must be taken into account. In addition care must 
be taken since the CMB probes the equation of state at 
a higher mean redshift and so may have a different effec- 
tive constant equation of state (eg. Saini et al. 2003) and 
is sensitive to perturbations. On the other hand there are 
some local measures that are largely insensitive to w or 
its evolution, such as the baryon fraction in clusters and 
mass-to- light ratios (Bahcall et al. 2000J. A recent analy- 
sis of the baryon fraction in clusters using Chandra data 
{Allen, Schmidt fc Fabian 2002^ puts a ten per cent uncer- 
tainty on the matter density when HST key project and 
nucleosynthesis information is used. 



4.4 Effect of priors on w 

The evidence is affected by the assumed choice of priors on 
the parameter values. Qualitatively, increasing the width of 
the priors on w will increasingly disfavour models with more 
parameters, which is clearly seen in Eqn.HO The priors in the 
previous section were arbitrarily chosen to be consistent with 
the input models and to have enough room for a fair degree 
of uncertainty prior to the data. We also note that although 
in principle one could imagine physical reasons for limiting 
Wo within a certain range, the dimensional quantity Wi = 
dw/dz could in principle be arbitrarily large. Therefore, a 
better, physically motivated prior is not on the polynomial 
coefficients but on the w(z) itself. For a subclass of dark 
energy models the equation of state satisfies the constraint 
— 1 < w(z) < 1. Figure |S] shows the effect of choosing this 
prior. We find that the evidence is better able to pick the 
correct model in this case relative to the choice of priors in 
the last section. For this subclass of theoretical models the 
prospects for disentangling the dark energy properties are 
seen to be much better. 
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Figure 5. Log of evidence plotted for the input model wq = —0.8 
and uii = 0.6 with f2 m = 0.3 as a function of f2 m assumed for 
the fit. The solid line corresponds to the fit to a cosmological 
constant, the dashed line corresponds to the fit to a constant w 
and the dot-dashed line corresponds to the linear fit respectively. 
It should be noted that the evidence is left un-normalizcd. 




Figure 6. Evidence values for various models, and various levels 
of fit as indicated on the x-axis, with a Gaussian prior fim — 
0.3 ± 0.05, and a uniform prior on the equation of state in the 
range — 1 < w(z) < 1. The dotted line is for an input model with 
w(z) = —0.8 + 0.6z, the short dashed line for an input model 
with w(z) = —0.7, the long-dashed line for a model with w(z) = 
-0.8 + 0.3Z. 



5 CONCLUSIONS 

We have investigated the use of Bayesian evidence for estab- 
lishing the significance of fitting cosmological data to various 
models of the dark energy. By focusing on a polynomial ap- 
proximation for the unknown equation of state parameter 
w(z) we have shown how evidence can be used to fix the 
polynomial order, thereby establishing the minimum varia- 
tion in w(z) required to fit the data adequately. 

We find that the evidence is affected by the priors. The 



largest uncertainty in establishing the nature of dark energy 
stems from the lack of knowledge of the precise value of the 
present day density in the form of pressureless dark mat- 
ter. Our results for the currently available SNe show that 
the current data favours the simplest case of cosmological 
constant. Since the peak of the likelihood lies in the range 
wo < —1, the tight prior —1 < w(z) < 1 disfavours the more 
complex models more strongly. 

We have shown that for a SNAP like data set, if one 
uses "wide" priors on the polynomial coefficients, evidence 
will enable us to decide if the data is best fit to a model with 
a cosmological constant, by a constant equation of state or, 
if the evolution is larger than wi ~ 0.5, to a linear model. 
We have shown that our conclusions can be significantly im- 
proved if we narrow the range of dark energy parameters by 
employing a prior of the form —1 < w(z) < 1, which hap- 
pens naturally for a subclass of dark energy models. This 
prior establishes the correct order of polynomial even for rel- 
atively modest evolution in the equation of state. However, 
since there are models which allow, albeit in a contrived way, 
w < —1, we have allowed for these models in our examples. 

To summarize, we have shown that evidence can be a 
powerful tool for systematically pinning down the nature of 
dark energy. Its application to the polynomial approxima- 
tion for the equation of state parameter allows us to fix the 
polynomial order that is required by data. The simple rela- 
tion between the polynomial coefficients and the underlying 
true equation of state enables us to obtain all the properties 
of dark energy provided by the data. 
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